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ABSTRACT 

A simple, observationally-motivated model is presented for understanding how halo masses, galaxy stellar 
masses, and star formation rates are related, and how these relations evolve with time. The relation between 
halo mass and galaxy stellar mass is determined by matching the observed spatial abundance of galaxies to 
the expected spatial abundance of halos at multiple epochs — i.e. more massive galaxies are assigned to more 
massive halos at each epoch. This "abundance matching" technique has been shown previously to reproduce 
the observed luminosity- and scale-dependence of galaxy clustering over a range of epochs. Halos at different 
epochs are connected by halo mass accretion histories estimated from A^-body simulations. The halo-galaxy 
connection at fixed epochs in conjunction with the connection between halos across time provides a connection 
between observed galaxies across time. With approximations for the impact of merging and accretion on the 
growth of galaxies, one can then directly infer the star formation histories of galaxies as a function of stellar 
and halo mass. This model is tuned to match both the observed evolution of the stellar mass function and 
the normalization of the observed star formation rate - stellar mass relation to z ~ 1 . The data demands, for 
example, that the star formation rate density is dominated by galaxies withM star ~ 10 M Q fromO < z < 1, 
and that such galaxies over these epochs reside in halos with M vlr 10 2 M . The star formation rate - 
halo mass relation is approximately Gaussian over the range < z < 1 with a mildly evolving mean and 
normalization. This model is then used to shed light on a number of issues, including 1) a clarification of 
"downsizing", 2) the lack of a sharp characteristic halo mass at which star formation is truncated, and 3) the 
dominance of star formation over merging to the stellar build-up of galaxies with M stal < 10 n M Q at z < 1. 
Subject headings: cosmology: theory — dark matter — galaxies: halos — galaxies: formation — large-scale 
structure of universe 



1. INTRODUCTION 

A fundamental goal of galaxy formation studies is to un- 
derstand what processes govern the stellar content and star 
formation histories of galaxies. A key piece of this puz- 
zle is relating the stellar masses and star formation rates 
of galaxies to the masses and formation histories of their 
associated dark matter halos. Ideally, one would like to 
make this connection by understanding the physical mech- 
anisms responsible for it from first principles. However, 
even the best current physically-motivated models of galaxy 
formation rely on significant approximations of unresolved 
physics. These approaches, based e ither on semi-analytic 
mode li ng (e.g. [White & FrenkHl99U I Somerville & Pr imack 
19991: IColeetalJ l20QOt lHatton etalJ l2003t ISpringel etaLl 
20011: ICroton et all l2006t Hiower et al.1 120061). or on hydro- 
dynamical simulations (e.g.jCe n & Ostriker 1992i;| Katz et alj 
ll996HSpringel & Hernquistl2003l:lKeres et alj2005l) still have 
trouble reproducing many basic observational results and suf- 
fer from serious uncertainties in the physical ingredients of 
the models. Although substantial progress has been made in 
these modeling efforts in recent years, star formation histories 
in these models and simulations are still sensitive to the inter- 
actions between a number of relatively unconstrained physical 
processes. 

Recent observation s have begun to measure the galaxy 
stellar mass functio n (Fontana et al. 20 04t iDrorv et al.1 1 2004 : 
Bundv et all 120051: iBorchet alJ 120061: iFontana et al.1 12OO6 ; 



Cimatti et al. 2006; Andreon 2006) and the star formation 



rate dNoeske et al.l2007btlZheng et al.l2007ah at high redshift, 
which complements more precise measurements locally (e.g . 
ICole et al.1 l200lt iBell et all l2003bt iBrinchmann et al.1 l2"004t 



Pante r et all l2007t ISalim et all 120071: ISchiminovich et alj 
120071) . At the same time, the evolution of dark matter halos, 
includ ing their abunda nce (e.g. | Warren et al.l2006l:lReed et al] 
1 20071). substructu res dKravtsov et al.1 12004 iGao et all 120041 
iReed et al.l 120051), a nd merger and accretion histories (e.g. 
Wech sler et al.l l2002). are becoming ever better understood in 
the context of the ACDM paradigm using numerical simula- 
tions. 

Several methods have recently been developed that take ad- 
vantage of these advances to connect the observed galaxy pop- 
ulation with dark matter halos using more empirical methods. 
The most popular of these, known as halo occupation mod- 
els, typically constrain the statistics of how galaxies populate 
their host hal os using galaxy clustering statistics and space 
densi t ies (e.g. | Scoccimarro et alj|2001l: iBerlind & Weinberg! 
I2002t iBullock et alJl2002t IZehavi et al.ll2004l) . An emerg- 
ing alternative is to connect galaxies to the underlying dark 
matter structure directly, under the assumption that the stel- 
lar masses or luminosities of the galaxies are tightly con- 
nected to the masses or circular velocities of dark matter ha- 
los. Throughout, this latter approach will be referred to as 
halo "abundance matching" because galaxies of a given stel- 
lar mass are matched to halos (including subhalos, which are 
halos that orbit within larger halos) of the same number den- 
sity or abundance. This approach matches the observed stel- 
lar mass function by construction, but has no other observa- 
tional inputs. Such an approach provides an excellent match 
to a number of galaxy clustering statistics at multiple epochs 



to a number ot galaxy clustering statistics at 
dKravtsov et al.N2004j: TTasitsiomi et al.1 2004 
1 20041: IConrov et al.ll2006l : lBerrier et alb oOf 
l2006llMarfn et al.ll2008l) . 
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The idea of abundance matching galaxies with dark mat- 
ter halos is not new, and it has been applied to associate a 
variety of objects wi th halos since the development of the 
CDM paradigm (e g. iMo et alj|1996t iMo & Fukugital 119961; 
ISteidel et ailll998t PWechsler et al.lll998l) . However, its suc- 
cessful implementation as a predictive tool requires a full ac- 
counting of the halo population, including the substructures 
that host galaxies, as well as a full accounting of the evolution 
of the abundance of galaxies as a function of their properties. 
These elements have only been in place quite recently. 

Halo occupation models as well as abundance matching 
models have been used primarily to understand the connec- 
tion between galaxies and halos at a fixed epoch, but recent 
work has begun to use these models to investigate the evolu- 
tionary history of galaxies, by combining information about 
the galaxy-halo connection at given epoc hs with theoretical 
input on the evolution of dark matter halos dWhite et aLll2007b 
IConrov et aUl2007bt IZheng et al.ll2007bt IConrov et alJl200a . 
In this paper we take the basic idea of abundance matching 
further, and use it to understand the evolution of the stellar 
content of galaxies. We use a simple, analytic representation 
of this framework, which connects dark matter halos to galax- 
ies by matching their abundances, to understand the build-up 
of stellar mass and the implied star formation rate of galaxies 
as a function of mass. We focus primarily on redshifts less 
than one, where the observational results are most reliable, 
but we expect the approach can be applied more widely and 
to earlier epochs as observational results improve. 

A complementary app roach has recently been presented by 
iDrory & Alvarez (2008). While we use the measured galaxy 
stellar mass function to connect galaxies to dark matter ha- 
los and infer the stellar mass buildup and star formation rates 
of galaxies, they used the measured star formation rates as a 
function of stellar mass, along with the time derivative of the 
galaxy stellar mass function, to infer the galaxy merger rate. 

The elements of our model are described in detail in ^21 32] 
presents our primary results, including comparisons to obser- 
vations. We discuss some of the implications of our model in 
$3] and summarize in Sj5] Throughout a flat, ACDM cosmol- 
ogy is assumed with the following parameters: (Cl m , 17a, a%) = 
(0.24,0.76,0.76), and h = 0.7 where h is the Hubble param- 
eter in units of 100 km s" 1 Mpc -1 . These cosmological pa- 
rameters are consis tent with the 3rd y ear WMAP estimates 
(Sper gel et ani2007h . A lChabrierl(l2003l) initial mass function 
(IMF) is adopted throughout. 

2. THE MODEL 

This section describes the details of our model. We start 
with a brief overview, and then move to a discussion of the 
halo ma ss function and galaxy stellar mass functions in £ 12.21 
and !2.3l The method used to assign galaxies to halos is out- 
lined in m.4\ followed by a description of the approach used 
to connect galaxies and halos across epochs in £12.51 Intro- 
ducing a sim ple e stimate for the effect of galaxy mergers and 
accretion in £|2.6l then allows us to compute star-formation his- 
tories of galaxies, as discussed in £12.71 

2.1. Overview 

The model described in detail in the following sections 
is an extension of previous modeling efforts that have been 
shown to succe ssfully reproduce an array of data from z ~ 
5 to the present (Kravtsov et al. 2004; Tasits iomi et alJ|2~004|; 



approach is to match the observed abundances of galaxies as a 
function of stellar mass with the expected abundance of dark 
matter halos. This step effectively assigns the most massive 
galaxies to the most massive halos monotonically and with 
no scatter. Since we include dark matter subhalos, which are 
halos orbiting within larger halos, we automatically include 
galaxies that would be observationally classified as satellites, 
although they are sub-dominant by number (~ 10-30% of 
the galaxies are satellites at any epoch). Thanks to parame- 
terizations of both the evolution of the observed galaxy stellar 
mass function and of the theoretical halo mass function, this 
connection between galaxies and dark matter halos can be de- 
termined continuously from z ~ 2 to z ~ 0. 

The novel feature of our approach, compared to previous 
work, is the use of average dark matter mass accretion histo- 
ries to connect the relations between halos and galaxies across 
time. jV-body simulations suggest that the aver age dark matter 
halo g rowth is a simple function of its mass dWechsler et al.l 
2002); thus, a halo at any given epoch can be connected to 
its typical descendants at later epochs. With the connec- 
tion between galaxies and halos determined at each epoch, 
the connection between halos across time implies an average 
connection between galaxies across time. At this stage the 
model produces the average stellar mass growth of galaxies 
as a function of both galaxy and halo mass. Since we use 
observationally-derived galaxy stellar mass functions as in- 
put, the connection is effectively one between observed galax- 
ies at different epochs. 

The final step is to differentiate these average stellar mass 
growth curves to infer the average mass-growth rates of galax- 
ies. The complication here is separating the growth due to star 
formation from that due to merging/accretion of other stel- 
lar systems. We introduce simple estimates of the contribu- 
tion due to merging that should bracket the possible effects 
of merging. This model then allows us to determine the av- 
erage star formation rates of galaxies as a function of their 
halo mass and redshift, which provides a key constraint on 
galaxy formation models. The following sections describe 
this framework in further detail. 

2.2. The halo mass function 

We use the cosmolo gy- and re d shift-d ependent halo 
mass function given by Warr en et al.l (120061) and transform 
their masses to M v ; r using an NFW (iNavarroet al] 119971) 
density profile with the concentration-mass relation from 
iBullock et al.l (1200 ll). assuming the updated model parameters 
given bv lWechsler et al.l (|2006). Our definition of the virial ra- 
dius corresponds to region with density contrast Ayk = 187r 2 + 
82jc-39x 2 w ith respect to the mean matter density, where 
x = Q(z)-1 (lBryan & Normanlll998l) . At z = 0, A vir = 337, 
and at high redshift A v h asymptotes to 180; 

The halo mass function provided by War ren et al.l (2006) 
only considers distinct halos, not the substructure within these 
distinct halos. Substructure as defined herein consists of ha- 
los whose centers are within the virial radii of larger halos, 
denoted subhalos. Distinct halos, in contrast, are those halos 
whose centers are not within any larger halos. 

We assume that the subhalo fraction is described by 



"sub 
"tot 



= 0.2 



0.1 



■z, 



(1) 



Vale & Ostriker 2004; Conrov et al. 2006; Berrier et al. 2006; 
Vale & Ostrikerfl2006[ |Marfn et al.ll2008l The first step in our 



independent of distinct halo mass, which provides a rea- 
sonable fit to data from simulations (see e.g., Fig. 1 of 
IConrov et all [2006 we don't include the moderate decrease 
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FIG. 1. — Evolution of the galaxy stellar mass function from z ~ 2 to 
z ~ 0. Our fiducial model for the evolution of the mass function (lines) is 
compar ed to the following observational resu lts from the literature: Bell et al. 
1 2003b, SDS S: circles). TPantei -et al] (2007, SDS S: diamonds) IDrory et alj 
120051 FDF). IBorch et al] )200a CO MBO- 1 7), IPerez-Gonzalez et al] (2001 
Spitzer), and Fontana et al] J2006L MUSI C). The disagreement between 
model and data at z = 2 is discussed in £|3.5I 



of /sub with increasing mass indicated by simulation data, but 
this would have a small effect on our results). Note that the 
subhalo fraction is defined with respect to the mass of the 
subhalos at the epoch of their accretion. This mass, rather 
than the present subhalo mass, has been shown to better cor- 
relate with observ ed galaxy properties (IConrov et alj l2006t 
iBerrier et al.ll2.006b . We thus derive an approximate halo mass 
function that includes both distinct halos and subhalos using 
this fraction. The results presented below are fairly insensitive 
to this fraction because it is small; we include it for complete- 
ness. Throughout, we refer to both distinct halos and subhalos 
as halos. 

2.3. The galaxy stellar mass function 

At each redshift, the number density (j)(M,z) of galax- 
ies with stellar mass M* is assumed to be described by a 
Schechter function, 



4>(M,z) = (f>*(z) 



M 

M*(z) 



exp 



M 

M*(z) 



(2) 



where the free parameters 4>*(z), a*(z) and M*(z) are, in prin- 
ciple, functions of redshift. We take the evolution of M*(z) to 
be: 



log[M*(z)/M ] = 10.95 + 0.17z-0.07z z 



(3) 



which is similar to the form advocated by iFontana et al.l 
(2006). Note that the evolution in M* implied from the above 
formula is mild at z < 2. Since the constraints on a* are weak 
at higher redshift, we assume for simplicity that it does not 
evolve: 

a* =-1.25, (4) 

which is consistent within the errors with available data to 
z ~ 2 ( IFontana et al.ll2006l) . 

The evolution of <fi*(z) raises a subtle but important issue. 
Various authors have measured 6* in redshift bins to z ~ 2 



and then proceeded to fit the observed <jf(z) to a function of 
the form 4>*(z) oc (1 + z) _/3 . However, it is clear that, modulo 
small evolution in M* and a* (and corrections due to stellar 
mass loss; cf. m.li , the time-derivative of (j)*(z) is simply 
the cosmic star formation rate (SFR) density. The functional 
form above, for typical adopted values of (3 = 1 - 3, results in 
an i ncreasing SFR d ensity at late times. This is not observed 
(e.g. lHopkinsll2004T) . 

In order to alleviate this tension, we have chosen to con- 
strain the evolution in <fi*(z) by requiring both that it repro- 
duce the observed evolution in the stellar mass function and 
that its derivative match the normalization of the observed star 
formation rate - stellar mass relation to z ~ 1 . After some ex- 
perimentation with different functional forms, we adopt the 
following evolution of (j>* : 



*(z) = 2x 10" 3 e 



-0.5z~ 



Mpc 



(5) 



Note that our parameterization is by no means unique. We 
have simply attempted to match the observed cosmic star for- 
mation rate density implied by our model and the observed 
stellar mass functions by adjusting the form of 4>*(z) - This 
functional form is similar to that given in IWilkins et al.l (T2008) 



-3 



as the best-fit 



who proposed <j>*(z) = 2.5 x 10" 3 e"°' 68z ' Mpc" 
to a variety of stellar mass function data. 

The zero-points of the Schechter parameters approximately 
reproduce the local set of parameters determined by a vari- 
ety of authors (ICo e et alJ2001t iBell et alJl2003bl: IWang et al l 



2006; iPanteret all 120071) . Figure [Tj compares the evolution 
of the galaxy stellar mass function in our model to various 
observational estimates. Our adopted Schechter parameters 
somewhat overpredict the abundance of low-mass galaxies at 
z = 0.5 and underpredict the abundances of all galaxies at z = 2. 
The latter disagreement is discussed in 0.51 

There is a second, perhaps more serious, tension created 
by comparison of the evolution of the stellar mass func- 
tion and the global SFR density. At z > 1 the integral 
of the star form ation rate density does not equal the stel- 
lar mass density (iNagamine et al] 20061: iHopkins & Beacoml 
I2006t IPerez-Gonzalez et al.l 120081; lwilkins et alj|2008l) . This 
tensi on can largely be removed if th e IMF evolves with red- 
shift dDavel2008l:IWilkins et alJ2008l) because the SFR probes 
the high-mass end of the IMF while the bulk of the stellar 
mass is contained in low-mass stars. An evolving IMF at z > 1 
is also suggested by recent work on the evolution of th e mass- 
to-light ratio of elliptical galaxies (Ivan D okkum 2008) a nd the 
abundance pattern s of metal-poor stars (Lucatello et al.1 l2005b 
lTurrdmsonll2007allbh . Whether this is the ultimate solution, or 
whether the solution lies in a more mundane systematic er- 
ror in one of the measured quantities is not currently clear. 
Because of this tension at high redshift, we focus our anal- 
ysis below z ~ 1, where the cosmic SFR density and stellar 
mass density are consistent with each other assuming a non- 
evolving IMF. 

2.4. Abundance matching: from halos to galaxies 

We assume that every galaxy resides in a dark matter halo 
and that there is a tight connection between the stellar mass 
of a galaxy and the mass of its associated dark matter halo. In 
the limit of zero scatter between galaxy and halo mass, halos 
of a given mass can be connected to galaxies of a given stellar 
mass by matching their abundances directly: 



n g (>M s ^ i ) = n h (>M vivJ ) 1 



(6) 
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i et alJ 
et alj 



where n g and «/, are the galaxy and halo mass functions, re- 
spectively (note that the halo mass function here includes both 
distinct halos and subhalos). In effect, this prescription as- 
signs the most massive galaxies to the most massive halos 
monotonically. Although the assumption of zero scatter is ide- 
alized, several recent works indicate that this scatter is small, 
with ~ 0.15 dex of s catter in galaxy luminosity at fixed mass 
( Zh eng et ai]|2007ri Ivan den Bosch et al.ll2Q07t lHansen 
120071; IWechsler et all I2008IT As shown in Tasitsiomi > 
(2004), scatter only effects the halo-stellar mass relation at 
the high mass end, and is in the sense that, at fixed galaxy 
mass, the mean halo mass decreases with increasing scatter. 

Since more massive halos are more strongly clustered at 
all epochs, this mapping implies that more massive/more lu- 
minous galaxies will also be more strongly clustered than 
less massive/less luminous ones, in qualitative agreement 
with a variety of clustering measurements ("e.g. IZehavi et alj 
12001 ICoiletail 12001 iLi et al.1 l2006t iMeneuxet all 120081) . 
This simple approach is surprisingly successful at quan- 
titatively matching an array of observations at multiple 
epochs and scales including mass-to-light ratios, cluster- 
ing measurements, and close pair counts dKravtsov et alj 
2001; iTasitsiomi et alj I200H IVale & Ostrikeri 120041 l2006t 
Berrier et all 120061; IConrov et all 120061; iMarin et all 120081) . 



confirming that it can be used with confidence herein. 

2.5. Mass accretion histories: from halo growth to galaxy 

growth 

Analysis of cosmological A^-body simulations has shown 
that ha lo mass growth can be described by a simple functional 
form dWechsler et al.li2002l) : 



M vir (a) = M D exp 



-2a, 



(t-0 



(7) 



where a = (1 +z)~ l , M is the mass of the halo at the redshift 
of observation a , and a c is the average formation scale factor 
of the halo, the single free parameter in the func tional form 
defined above. Following Wechsl er et al.l (|2002), we adopt 
the following parameterization of a c , which provides a good 
fit to jV-body simulations: 



fl c (Mv5r) = 



4.1 



c(M vir )(l+z) 



(8) 



where c(M v ; r ) is the halo concentrat ion- halo mas s rela tion at 
z = 0. We use the model given by Bullock et al. (200B and 
the updated parameters provided in Wechsler et all (120061) for 
c(M vir ). 

In the previous section we showed how to construct M star - 
M v ir relations as a function of redshift. The relation between 
a halo of mass M v ; r at one epoch and its mass at some lat- 
ter epoch is known via Equation [7] This relation between 
halos across time allows us to connect the M star -M v i r rela- 
tions across time and hence allows us to determine the stellar 
growth of galaxies. 

For example, we can start with a halo mass M V u at some 
early epoch. The M stal -M wa relation at that epoch then deter- 
mines the stellar content of the halo. We can then evolve this 
halo to a later epoch via Equation|7] With the A/ Star — rela- 
tion at this later epoch we can then read off the stellar content 
of the halo at this later epoch. Continuing this process allows 
us to build up the full stellar mass growth of the galaxy sit- 
ting at the center of this evolving halo. This process can be 
repeated for all halos of all masses, allowing one to determine 



the stellar mass growth of galaxies as a function of dark matter 
halo mass. 

Equation [7] does not apply to subhalos and yet we have in- 
cluded subhalos in our analysis up to this point. There are at 
least two reasons why this issue will not significantly impact 
our results. First, at any given epoch the major i ty of subhalos 
were only recently accreted dGao et al.1 120041; IZentner et al.1 
120051) . and thus Equation [7] should provide a reasonable ap- 
proximation to the mass growth history of subhalos over 
most of their evolution. Second, as mentioned above, sub- 
halos constitute a small fraction of the total halo population 
(~ 10-30% at any epoch) and thus this approximate treat- 
ment should have a small effect on our conclusions. 

2.6. The impact of merging on galactic growth 

Galaxies can grow in stellar mass by either star for- 
mation or by the cannib alism of smaller galaxies (e.g. 
lOstriker & Hausrrianlll977l) . Both processes can in general 
contribute to the average stellar mass growth of galaxies. We 
are interested primarily in the inferred star formation rates as 
a function of redshift, stellar mass, and halo mass, and we 
thus seek a simple way of accounting for the impact of galaxy 
merging on the stellar mass growth of galaxies. In what fol- 
lows we present an estimate for the accretion rate of smaller 
galaxies onto the halos of larger galaxies. We then consider 
two assumptions for the fates of these accreted systems that 
will bracket the range of possibilities. In the first, we allow 
all of the stellar material accreted onto the halo to rapidly 
merge with the central galaxy, thereby increasing the mass 
of the central galaxy. The other possibility we consider is 
that the accreted material remains within the host halo of the 
central galaxy but does not add to its measured luminosity. 
In other words, the accreted material either remains as bound 
satellite galaxies or ends up as diffuse stellar material not de- 
tected in standard survey photometry. In this latter scenario 
stellar mass growth is thus determined entirely by star forma- 
tion. These two scenarios are referred to as the "merger" and 
"no-merger" scenarios below. We now describe the merger 
scenario in more detail. 

Halos grow via the accretion of smaller halos. The mass 
spectrum of accreted halos is approximately self-similar in 
m' = mjM z where m and M z are the mas s of the accreted halo 
and mass of the par ent halo at redshift z dLacev & Coldfl993l; 
IStewart et al.ll2008l) . The spectrum can be approximated as: 



df 
dlnm' 



2.6 



exp 



m 
0J 



(9) 



where / is the fraction of mass accreted in clumps of mass ml . 
The exponential cut-off is steep because m 1 > 1 is not allowed 
by definition. Equa tion [9] is a fit to the simulation results of 
Stewa rtet al.l d2008). This function does not integrate to unity, 
indicating that a significant fraction, ~ 30-50%, of the par- 
ent mass is accreted in a diffuse component of dark matter 
(IStewart et a l. 2008). Whether or not this compo nent is truly 
diffus e or is in clumps of very small mass (e.g. lMadau et al.1 
2008) is immaterial for our purposes, because in either case 
this component will not bring in additional stars. 

With the mass accretion spectrum in hand, the halo growth 
rate, m^o, can be converted into a stellar growth rate due to 
mergers, m stals , via: 



^stars — '^halo 



df Mst 
dlnm' M v 



-(M vin z)dW, (10) 
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where jf-iM^i^z) is the redshift- and halo mass-dependent 
stellar-to-halo mass ratio determined in ^2.41 EquationfTOlcan 
be thought of as a convolution of the halo mass accretion spec- 
trum with the relations between halo mass and stellar mass 
determined in previous sections. In other words, Equation [9] 
tells us the types of halos that are accreted, and ^2.41 tells us 
the stellar content of these accreted halos. 

In the following sections we will present results for both the 
merger and no-merger scenarios — two scenarios that bracket 
reality. We remind the reader that while we discuss these sce- 
narios separately, in reality there is clear evidence that both 
cases occur. In particular, there exists direct observational ev- 
idence for galaxy mergers, and the existence of gravitation- 
ally bound groups and clusters of galaxies implies that not 
all galaxies merge when their halos merge. Thus, when in 
later sections we discuss a preference for one scenario over 
another, we do not mean to suggest that the other scenario 
never occurs but rather that it is of sub-dominant importance 
when attempting to describe the statistical properties of galax- 
ies. 

In i j4.3l we discuss how our results bear on the relative im- 
portance of merging and star formation on the stellar growth 
of galaxies. 

2.7. From galactic growth to star-formation rates 

With an estimate for the amount of stellar mass growth that 
is due to merging/accretion, the SFR of an average galaxy can 
then be estimated straightforwardly via a derivative of the por- 
tion of stellar mass growth attributed to star formation. The 
relation between mass growth and star formation is compli- 
cated by mass loss due to dying stars. We take into account 
this effect with the following formula: 

, /f + 3xl0 5 \ 
/ loss (f) = 5xlO- 2 lnf 3xiQ5 j, (11) 

where t is in years and f\ oss (t ) is the fraction of mass lost by 
time t for a co-eval set of stars. This formu la is a fit to the 
mass loss of simple stellar populations with alChabrierj (2003) 
IMF dRenzini & Ciottilll993UBruzual & CharlotH2003l) . Note 
that only ~ 60% of the stellar mass formed in a burst of star 
formation remains after several gigayears, and that the stellar 
mass remaining includes stellar remnants. With the full stel- 
lar mass growth curve one can then iteratively solve for the 
star formation rate required to generate such growth given the 
above mass-loss rate formula. 

As mentioned in ^2.31 the form we have chosen for the red- 
shift evolution of the stellar mass function is not unique. Yet 
it is clearly the rate of evolution of the mass function that de- 
termines the resulting star formation rates of galaxies in our 
model. We re-emphasize that this form was chosen to best 
match the observed cosmic star formation rate density. Our 
goal is not to find a unique form for the evolution of the mass 
function but rather to present a consistent framework in which 
to interpret a vast array of observational data, and to link that 
data to the underlying dark matter skeleton. 

3. MODEL IMPLICATIONS 

The previous section presented a method for connecting the 
stellar masses and star formation rates of galaxies to dark mat- 
ter halos over a range of epochs. This section explores these 
relations and compares to observations where possible. 

3.1. Halo- galaxy connections 



Implementing the abundance matching technique discussed 
in £|2.4l at various epochs yields the relations between halo and 
galaxy mass shown in the top panel of Figure [2] The generic 
shape of the relation is governed simply by the Schechter-like 
functions of both the galaxy stellar MF and the halo MF. The 
redshift outputs are spaced equally in (1 +z)~ l - One novel con- 
clusion drawn from this figure is that, since z ~ 2, the stellar 
mass of galaxies residing in halos of mass ~ 10 125 M Q stays 
roughly constant, at ~ 10 n M Q . Over the redshift range con- 
sidered, above this mass scale, halo growth out-paces stellar 
mass growth, while below this scale, galaxy growth is more 
vigorous than halo growth. 

The relation shown in Figure |2] applies to central galaxies 
and to satellites where the halo mass refers t o the mass at 
the ti me of accretion onto their host (see, e.g. IConrov et al.l 
I2006L for a discussion). Including a modest l evel of scatter 
between stellar and halo mass (as discussed in $2.41 does not 
substantially impact the mean relation shown in this plot. As 
a comparative reference, we show the location of the Milky 
Way in t his figure, as determ ined from the halo mass esti- 
mates of lKlvpinetal.1 d2002l) . The Milky Way falls directly 
on our mean relation. 

The shape and mild evolution of the M stal -M v i r relation 
shown in Figure [2] provides a clear interpretation of the ob- 
served relation b etween stellar and hal o mass from z ~ 1 to 
z ~ reported in IConrov et al.1 (l2007ch . These authors used 
the dynamics of satellite galaxies orbiting around brighter 
host galaxies to constrain halo masses, and found that in bins 
of galaxy stellar mass, halo mass evolves little or not at all 
since z ~ 1 below M star < lO n M but increases by a factor of 
several above this stellar mass. This qualitative trend is evi- 
dent in Figure |2j and can be attributed to the fact that above 
M star ~ 10 Mq the relation shallows, implying that a small 
shift in the relation over time produces a large change in the 
halo mass of a given galaxy mass over time. 

Usin g a similar approach to our abundance matching tech- 
nique, Sh ankar et alj (120061) find comparable results on the 
evolution of the M star - M v ; r relation since z ~ 1 . More- 
over, group catalogs constructed from large observational sur- 
veys have begun to probe the M star -M v j r relation at z ~ 
dBerlind & others.1120061: lYang et alj|2007l) . Results from the 
catalogs are in good agreement with what we find at z ~ 0. 

We define the integrated efficiency of star formation as 
rj = Afsta r/Myir/./fc, where fb = 0.17 is the universal baryon 
fraction dSpergel et al.1 120071) . This efficiency quantifies the 
fraction of available baryons that have been converted into 
stars, and peaks where integrated star formation is most ef- 
ficient. Abundance matching readily predicts r)(M wu ) and is 
shown as a function of redshift in the bottom panel of Fig- 
ure [2] The first thing to note is that this result implies that 
the overall efficiency of converting baryons into stars is quite 
low, never reaching more than ~ 20% of the potentially avail- 
able baryons. Although perhaps somewhat surprising, note 
also that the global stellar mass density is ~4-8 times less 
than the global baryon density. This low efficiency is also in 
good agreement with curr ent estimates for the total and stellar 
mass of the Milky Way (Kly pin et al.ll2002T) . with estimates 
of halo masses from weak l ensing measurements com bined 
with stellar mass estimates dMandelba um et ail 20061). with 
halo oc cupation models both at z ~ and z ~ 1 dZheng et al.l 
2007b), and with esti mates from the accounting of baryons 
in various states (e .g. iFukugitaet 111119981; iBell et alj|2003at 
iBaldrv et al.ll2008l) . 

Note that this quantity is not the instantaneous star forma- 
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FIG. 2. — Top Panel: The relation between galaxy stellar mass and halo 
mass from z = 2 to z = 0, using the abundance matching model. Bottom Panel: 
Fraction of available baryons that have turned into stars (integrated star for- 
mation efficiency) as a function of the halo mass and redshift, where ft is the 
universal baryon fraction. The star marks the location of the Milky Way at 
2 = 0. The thick black line represents the relation at z = 1. The relations at 
z > 1 (dashed lines) should be treated with caution; see 93.5l for details. 
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FIG. 3. — Top Panel: The redshift-dependent relation between galaxy stellar 
mass and the mass of the halo to which it will belong at z, = 0. By focusing 
on a fixed z = halo mass, one can read up in the plot from z = 2 to z = to 
infer the stellar mass growth of an average galaxy that, by z = 0, resides in 
that halo. Bottom Panel: Fraction of available baryons that have turned into 
stars as a function of the z = halo mass. The thick black line represents the 
relation at z =1. The relations at z > 1 (dashed lines) should be treated with 
caution; see Sj33]for details. 



tion efficiency because 77 herein is defined with respect to the 
stars that still exist in the galaxy. For the IMF we adopt, ap- 
proximately 40% of the stellar mass that forms is rapidly lost 
as massive stars die. One may convert 77 into the fraction of 
baryons that have ever spent time in a star by dividing the 
numbers we quote by 0.6. Thus, for example, the peak shown 
in the bottom panel of Figure [2] would be ~ 30% if one were 
interested in the instantaneous star formation efficiency. 

Two important trends are apparent. First, the location of 
the peak decreases to lower halo masses with time. The latter 
trend i s a manifestation of at least one meaning of "downsiz- 
ing" dCowie et al.lll996l) . and is a natural implication of the 
fact that the characteristic stellar mass evolves more slowly 
than the characteristic halo mass (see £14.1 1 for a discussion 
of downsizing). Second, the amplitude of the peak star- 
formation efficiency increases with decreasing redshift, al- 
though the magnitude of this trend is somewhat uncertain. 
This is due to the rather uncertain evolution of <fi*, which di- 
rectly affects the evolution of the peak of rj(M wu ); varying the 
evolution in </>* over a reasonable range changes the amount 
of evolution in the peak by less than a factor of two. More ac- 
curate observational constraints on <fi*(z) are required to more 
robustly pin down the evolution in 77(M v i r ). However, the trend 
that the peak shifts to lower masses with time is robust to un- 
certainties in </>*(z). The mass at which baryons are most ef- 
ficiently converted into stars shifts by about a factor of ~ 20 
fromz = 2toz = 0, fromM vir ~ 10 13 M Q toM vir - 1O 1L7 M . 

These trends, including the factor of ~ 2 increase in peak 
efficiency from z ~ 1 to z ~ 0, and the mild decrease in halo 
mass at which the peak occurs, agree well with halo oc- 



cupation modeling of galaxies at these epochs dZheng et alj 
2007bh . 

As discussed in ^2.51 we know from A^-body simulations 
the full mass growth histories of dark matter halos statistically 
for a given cosmology. These accretion histories allow us to 
evolve a halo of mass M v ; r (z) at redshift z forward in time 
to the mass such a halo will have by z = 0, M v - n -(z = 0). We 
can couple these halo accretion histories to the M stal -M v j, 
relations discussed above to determine the relation between a 
halo's mass at z = and the stellar mass content of that halo 
as a function of redshift. 

The resulting relations are shown in the top panel of Fig- 
ure [3] The z = halo mass can be thought of as a unique tag 
for each (average) galaxy. A vertical slice through Figure [3] 
thus traces out the trajectory of an average galaxy at differ- 
ent redshifts. From this plot one can thus read off directly 
the average stellar mass build-up of galaxies as a function of 
their z = halo mass. These relations are uniquely determined 
by the relations between galaxies and halos at fixed epochs in 
conjunction with the evolution of halos demanded by our fidu- 
cial cosmology. 

A clear and robust inference from this figure is that 
the stellar mass of galaxies residing in z = halos of 
mass > 10 14 M Q was mostly assembled by z ~ 2. This 
agrees qualitatively with the modest evolution in the mas- 
sive end of the observed galaxy ste l lar mass function 
since z ~ 2 (e.g. iFontana et all 120041: iDrorv et al.1 1 2004 : 



Bundv et al. 



Cimat tiet al.l 



20051: iBorch et alj 120061: Fontana et all 2006 : 
20061: lAndreonll2006t iBrown et al.ll2007l 12008 : 
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FIG. 4. — Fraction of mass assembled as a function of redshift. Each panel 
compares the fraction of a galaxy's stellar mass that is assembled (solid and 
dotted lines) to the fraction of the galaxy's parent halo mass that is assembled 
(dashed lines). The stellar growth curves are dott ed at z > 1 to indicate that 
this regime should be treated with caution; see ij3.5l for details. The four 
panels display the assembly history for a range of galaxies with z = stellar 
and halo masses shown in the legend, in units of log(Mo). It is clear that, at 
lower mass, a larger fraction of the halo is in place at early times compared 
to higher mass. The opposite trend is true for stellar masses. 



FIG. 5. — Average formation times for galaxies (solid line and shaded re- 
gion) and halos (dashed and dotted lines) as a function of z = galaxy and 
halo mass. Lines and shaded region indicate the redshift at which 25%, 50%, 
and 75% of the final mass was assembled, as labeled in the figure. For exam- 
ple, a galaxy with stellar mass 1O 1O M at z = assembled 50% of its final 
mass by z ~ 0.5. Such a galaxy resides in a halo of mass 10 1L6 M©, which 
was half assembled by z ~ 1.0. At /Vf st ar ^ 10 10 ' 7 Mq halos are assembled 
before galaxies while at higher masses the opposite is true, in agreement with 
the trends seen in Figure [4] 



iPerez-Gonzalez et alJl2008t ICool et all feOOS). Note however 
that the input to our model is the observed galaxy stellar mass 
function, so for example if observations do not account for 
the low surface brightness intracluster light associated with 
central galaxies in massive halos (e.g . [Gonzalez et al1l2005b 
Zibett i et al1l2005l: iKrick & Bernsteinll2007l) . then our model 
will also fail to incorporate this component. The diffuse light 
could contain as much mass as the central galaxy, and should 
thus be taken int o account when modeling massive galaxies 
(see discussion inlMonaco et alj|2006t IConrov et al.ll2007blal : 
lPurcelletal.ll2007tr Here however our focus is on stellar 
growth and star formation in more modestly-sized halos. 

The bottom panel of Figure [3] shows the integrated star for- 
mation efficiency of galaxies at various epochs as a function 
of their z = halo mass. For halo masses less than ~ 10 12 M Q 
the integrated efficiency is a monotonically increasing func- 
tion of time. At higher masses the efficiency rises, peaks, 
and then falls with increasing time, and at masses greater than 
~ 10 14 M Q the integrated efficiency is a continually decreas- 
ing function of time since z = 2. 

3.1.1. Galaxy growth versus halo growth 

A perhaps more revealing illustration of the results in Fig- 
ure [3] are shown in Figure |4] There the stellar mass growth 
of a galaxy is compared to the growth of its parent dark mat- 
ter halo for four representative z = stellar masses. For stel- 
lar masses < 10 10 ' 7 M o (corresponding to M vil - < 1O 12 3 M ), 
fractional halo growth since z = 2 is much more mild than 
stellar growth. This low-mass regime can thus be thought 
of as 'internally-dominated', where growth is not controlled 
by extra-halo processes. Stellar growth in this regime is thus 
driven by gas physics related to cooling, star formation, and 
feedback. 

The situation is qualitatively different at higher masses. At 



stellar masses > 10 n M Q (corresponding to M vir > 10 13 M Q ), 
fractional halo growth is much stronger than stellar growth at 
z < 2. This is not surprising in light of the fact that the massive 
end of the stellar mas s function appears to be approximately 
in place since z ~ 1 dFontana et al.ll2004t iDrorv et al J [2.0041 : 
Bundv et al.1 120051: iBorch et alj l2006t iFontana et al l 200fj 
Cimatti et all 120061: IWake et alj [2006; Br own et alj 120071: 
Cool et alJ l2008). This high-mass regime is thus 'externally- 
dominated', in contrast to lower-mass systems. High-mass 
systems are primarily accreting copious amounts of dark mat- 
ter, some fraction of which will bring in bound stellar systems 
(e.g. satellite galaxies). The system — defined loosely as 
the region within the halo virial radius — is thus growing in 
a larger sense, while the galaxy at the center of the halo is 
not. This figure does not include the potentially massive com- 
ponent of stellar light associated with a diffuse background, 
known as the intracluster ligh t. For a discussion of the impor- 
tance of this component see IConrov et al.l d2007bl) and refer- 
ences therein. 

The average formation times for galaxies and halos is 
shown in Figure [5] Here we plot the redshift at which 25%, 
50%, and 75% of the final mass was assembled, for both the 
stars accreted onto the galaxy and for the dark matter mass 
accreted onto the halo. Since we are investigating only aver- 
age properties, recall that in our model each galaxy mass is 
assigned a unique halo mass, and thus each average galaxy 
has a unique halo mass and stellar mass, and unique forma- 
tion times for both of those components. It is clear that the 
stellar mass in lower mass systems formed later than in high 
mass systems, while the dark matter halos display the oppo- 
site trend. Although this general trend has now been evident 
from a range of data, this figure ties together the available in- 
formation on galaxy and halo growth. 
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FIG. 6. — Top Panel: Average specific star formation rate history as a func- 
tion of redshift, for galaxies with z = stellar masses given in the legend (in 
units of log(Mo)). The dotted line indicates a slope of unity. Bottom Panel: 
Star formation rate history as a function of redshift for the same model galax- 
ies. The solid and dashed lines represent our model for two prescriptions to 
relate stellar growth to star formation. The solid lines are f or the no-merger 
model while the dashed lines are for t he merger model (s ee i|2.6l for details). 
Observed star formation histories from Panter et al. (2007) for galaxies of the 
same z = stellar masses as the model predictions are included for compari- 
son. For reasons discussed in the text, the model is likely not reliable at z > 1 . 
Note that results for both the model and data are average relations. 



Notice that the results presented thus far do not require 
making any assumption for how the galaxy mass is built up 
(i.e. whether by star-formation or merging). These results 
only require knowledge of the redshift-dependent galaxy-halo 
connections in conjunction with how average halos grow with 
time. We now turn to the more challenging task of constrain- 
ing the possible modes of galaxy growth. 

3.2. The star formation history of galaxies 

The star formation rate (SFR) and specific star formation 
rate (SSFR = SFR/M star ) for average galaxy trajectories are 
shown in Figure [6] as a function of redshift for several rep- 
resentative z = stellar masses. The SSFR is approximately 
self-similar in galaxy mass for z < 1 and scales with redshift 
as log(SSFR) oc z. Lower mass galaxies have higher specific 
star formation rates at all times. In this and subsequent fig- 
ures we include two different treatments for the importance of 
merging on stellar mass growth. The solid lines represent the 
assumption that all stellar mass growth is due to star forma- 
tion in situ, while the dashed lines represent the prescription, 
described in ^2.61 that includes mergers and accretion. 

The SSFR for the most massive galaxies (M star > 10 n M o ) 
is the least constrained in our model because the massive end 
of the stellar mass function evolves little at z < 2, and so the 
SFR is a derivative of a nearly constant function. This can be 
seen in the top panel of Figure [3] where it is clear that small 
uncertainties/changes in the stellar mass or halo mass scale 



can cause relatively large uncertainties in the star formation 
rates. When discussing star-formation rates we restrict our- 
selves to stellar masses less than ~ 10 n M Q where the con- 
version between stellar mass growth and SFR is most reliable. 

The bottom panel of Figure [6] shows the average SFR of 
galaxies as a function of z = stellar mass. The results in this 
figure are not directly accessible to observations at high red- 
shift because the observations do not tell us the connection be- 
tween galaxies at different epochs. Without this information, 
following the SFH of a particular galaxy across time is not 
possible with observations of the distant past. However, such 
information is at least in principle attainable from consider- 
ation of the stellar populations of galaxies in the local uni- 
verse. Attempts at inferring the SFH of galaxies in this way 
have concluded that the SFR peaked earlier for more mas- 
sive galaxies, and for galaxies less massive than - 10 10 M o 
the dat a are consistent with a constant SFH, at least since 
z ~ 1 dLee et alj|2007t Ivan Zed 120011: iHeavens et all l2004t 
iPanter et al.ll2007l) . The model results presented i n Figure [6] 
is compared to the results from lPanter et alJ (120071) who have 
used the stellar populations of local galaxies to constrain their 
mass-dependent star formation histories. Our model repro- 
duces the general trends well, though there is a ~ 0.1-0.2 
dex offset between the model and data. 

Furthermore, the trends in Figure [6] sh ed light on the phe- 
nomenon known as 'downsizing' ("e.g. ICowie et al.l 119961 : 
iBrinchmann & E llis 2000: lJuneau et al.ll2005l) whereby more 
massive systems formed the bulk of their stars at earlier 
epochs compared to less massive systems. While it is clear 
that star formation has peaked earlier in more massive sys- 
tems, it is also apparent from the figure that at any epoch, 
more massive galaxies have higher star formation rates and 
smaller specific star formation rates than less massive galax- 
ies. We discuss this issue further in 34.11 



3.3. SFR dependence on galaxy and halo mass 

Our model allows us to calculate the SSFR and SFR of 
galaxies as functions of stellar mass and redshift; these are 
shown in Figure [7] As above, the figure includes two differ- 
ent prescriptions for relating stellar growth to star formation. 
One possibility is that all stellar growth is due to star forma- 
tion (solid lines) while the other allows for some fraction of 
stellar growth to be attributed to merging {dashed lines; see 

The model is compared to a variety of data from the lit- 
erature over the redshift interval < z < 1 . In all cases the 
data are meant to represent average star formation rates as 
a function of stellar mass (i.e. the average star formation 
rate of all — both red and blue — galaxies at a given stel- 
lar mass). T he completeness-corr ected average SFR-Mstar 
relation from lNoeske et al . ( 2007b) was construct ed based on 
the completeness corrections of iLin et a l. (2008, K. Noeske 
private communication) in order to account for red galaxies 
with no detectable l evels of star formation. The results from 
Zheng_etalJ (l2007al) were derived from stacked data and can 
thus b e interpreted as a verage relations. Finally, the results 
from lSalim et alJ d2007l) were determined from data with suf- 
ficient sensitivity to detect extremely low levels of star for- 
mation and can thus also be interpreted as an average relation 
over all galaxies at a given stellar mass. 

It is clear from Figure [7] that the assumption that all stellar 
growth is due to star formation (i.e. the no-merger scenario; 
solid lines) provides a much better match to the data at stel- 
lar masses > 10 10 M Q . For this reason we adopt this assump- 
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FIG. 7. — SSFR (top panel) and SFR (bottom panel) as a function of galaxy 
stellar mass at various epochs. Note that t he two pan e ls displ ay equivalent 
information. Observational results from Noeske et al. 1 2007b) are included 
and labeled N07, as are z ~ data from Salim et al. ( 2007) and data from 
Zheng et al. (2007a), labeled Z07. Error bars denote Icr scatter, not the error 
on the mean, and are included for only one set of data for clarity. The solid 
and dashe d lines represent the no-merger and merger models, respectively 
(see i]2.6l for details). The no-merger model is clearly a better match to the 
data at all epochs. Note that the model describes average relations of SFR 
with mass. 



tion as the fiducial model. At lower masses the no-merger and 
merger scenarios yield the same predictions for the star forma- 
tion rates (i.e., even with maximal merging, incoming satellite 
galaxies do not contribute any appreciable stellar mass) The 
importance of star formation over merging in galactic growth 
is discussed further in 34.31 

It is instructive to understand which aspects of the model 
are driving agreement with the data. The normalization of the 
model predictions depend on the evolution of <fi*, the redshift- 
dependent normalization of the stellar mass function, while 
the shapes depends on a* and M* , although the latter two 
dependencies are much weaker than the first. Recall that we 
have tuned the evolution of <fi* to reproduce the normalization 
of the SFR-M star relations, but not the shape of these relations. 
The shape is thus a robust prediction of our approach, while 
the normalization agrees with the data by construction. 

Figures [6] and [7] can be thought of as consistency checks 
between the model and data, since many of the implications 
that can be drawn from these figures are readily available from 
the data themselves. In contrast, Figure [8] contains a variety 
of novel results. This figure shows the model predictions for 
the SFR and SSFR in galaxies as a function of their host dark 
matter halo masses. 

An interesting consequence of Figure [8] is that the halo 
mass at which the most vigorous star formation occurs is 
not a strong function of redshift. In addition, the peak in 
SFR occurs over a large range of halo masses, rather than 
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FIG. 8. — SSFR (top panel) and SFR (bottom panel) as a function of halo 
mass at various epochs. As in previous figures, the solid and dashe d lines rep- 
resent the no-merger and merger models, respectively (see !j2.6l for details). 
These two assumptions only impact the derived SFR at high masses, and the 
data favor the no-merger model (see Figure[7}. In the top panel the dotted line 
indicates a slope of —0.5. Note that the model describes average relations of 
SFR with mass. 

at one well-defined scale. Moreover, the SSFR-M v j, rela- 
tion appears to be almost scale-free up to M v ; r ~ 10 13 Mq 
for our favored model, with the normalization steadily de- 
creasing with time and a non-evolving slope. Over the range 
10 110 <M vir < 10 13 M Q , the redshift- and mass-dependent 
relation can be approximated by: 

SSFR » 1() ; " Gyr" 1 . (12) 

The relation between SFR and halo mass shown in Figure[8] 
can be thought of as the most fundamental of the relations dis- 
cussed in this work, as this redshift-dependent relation gives 
rise to all other relations. This relation thus provides a di- 
rect link between observations and models in the sense that 
any model which reproduces the trends in Figure [8] will auto- 
matically match the variety of observational results discussed 
herein. This connection between SFR and halo mass, deter- 
mined entirely from observations with our simple approach, 
can thus be of general use to the modeling community in con- 
straining models of cooling, feedback, and star formation in 
galaxies. 

Figure [9] shows the SFR-M v h relations again, now with 
Gaussian fits (note that the y-axis is shown here in linear 
units). It is clear that the relations are well-characterized as 
Gaussian except perhaps at the highest masses at low red- 
shift (where our model is least well-constrained) and at high 
masses at high redshift where the turn-over is not resolved. 
The fits in these regimes should thus be treated with caution. 

The Gaussian fits are characterized by three parameters: the 
peak, SFRo, mean, Mq, and dispersion, a. The first two pa- 
rameters are shown as a function of redshift in the inset pan- 
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FIG. 9. — SFR as a function of halo mass at various epochs for our preferred, no-merger model {solid lines). In this figure we include Gaussian fits to the 
SFR(M v j,) relations {dashed lines). The inset shows the best-fit normalization, SFRo, and mean. Mo, as a function of redshift {solid lines), along with linear fits 
{dashed lines) to these relations, which are barely distinguishable from the relations themselves. The dispersion is a weak function of redshift and is thus only 
included in the three larger panels for clarity — the average dispersion between 0. 1 <z< 1 is a = 0.72. Note that at higher redshift the turn-over at high masses 
is not resolved and the fits there should thus be treated with caution. 



els of Figure[9] These parameters are very well approximated 
with the following linear relations: 

log(SFR ) = 0.47+ 1.1 z (13) 

log(M )= 12.3 + 0.81 z. (14) 

These fits are included in the inset panels. The dispersion is 
a much weaker function of redshift than either the normaliza- 
tion or the mean. The dispersion ranges from 0.64 at z = 0.1 
to 0.89 at z = 1.0 with a mean value of 0.72 over the whole 
interval 0.1 < z < 1.0. It is important to recognize that, while 
the general functional form and redshift-dependent trends are 
robust predictions of our model, the precise values are subject 
to uncertainty because the observations themselves, to which 
the model is tied, still have substantial uncertainties. 

In all of these figures it is important to keep in mind that we 
are presenting average relations between various quantities. 
At first glance Figure [8] might suggest that there would be no 
red galaxies (where star-formation has ceased) at z ~ 1 . There 
can of course be such galaxies, as there can also be galaxies 
with SFR in excess of the average relation presented in Figure 
i 

It is worth mentioning here why the approach taken in this 
paper is particularly useful. The current generation of hy- 
drodynamic simulations and semi-analytic models are not ca- 
pable of re producing the redshift-dependent trends shown in 
Figure □ (Dave 2008). The c ause of this discrepancy is not 
currently understood, although Dave (2008) speculates that an 
evolving IMF can alleviate the tension. Regardless, it is clear 
that until this tension is resolved, using either hydrodynamic 
simulations or semi-analytic models to interpret the observa- 
tions and connect them to the formation and evolution of ha- 



los requires caution. Our approach matches the observations 
by construction and it can thus be used with more confidence 
for interpreting the data. Its main limitation, and the main 
advantage of the simulations, is that our model makes no ref- 
erence to the underlying physical processes governing these 
relations. However, the connection between observables and 
halo mass derived from our approach should be very helpful 
in informing these more physical models. 

3.4. Global properties 

FigurefTUlplots the SFR density as a function of stellar mass 
at z = 0.1,0.5, and 1.0. This quantity is the SFR density con- 
tributed by galaxies with mass M stal - and is produced by multi- 
plying the SFR-M stal relation by the galaxy stellar mass func- 
tion, $(M sta r)- This quantity is thus well-constrained by the 
observational data. In the figure, this quantity is plotted both 
as a function of the stellar and dark matter halo mass. 

In the top panel, the peak for our favored model (solid lines) 
is ~ 0.5 dex lower than the characteristic mass of the stellar 
mass function, M*, at all epochs, indicating that the bulk of 
the SFR density is contributed by galaxies a factor of ~ 3 in 
mass below M*. In other words, the characteristic galactic 
mass in which stars form since z < 1 is a factor of ~ 3 lower 
than the characteristic galactic mass dominating the mass den- 
sity. 

This peak in the SFR density does not change appreciably 
from z ~ to z ~ 1. The bulk of star formation at z < 1 
never occurs in small systems, rather it is always dominated 
by relatively massive, M star - 1O 1O °- IO ' 5 M galaxies, simi- 
larly, the peak as a function of M vn - (bottom panel) decreases 
only slightly, by at most ~ 0.5 dex from z = 1 to z = 0. These 
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FIG. 10. — Top Panel: SFR density as a function of stellar mass and red- 
shift. Bottom Panel: SFR density as a function of halo mass and redshift. 
These plots illustrate the contribution to the global SFR density for galaxies 
of a given stellar mass (top panel) and for galaxies residing in a given halo 
mass (bottom panel). As in previous figures, we include SFR estimates for 
both the assumption that the stellar growth is entirely due to star formation 
(solid lines), and a simple prescription to account for the amo unt of stellar 
growth due to mergers and accretion (dashed lines, see i]2.6t . This figure 
demonstrates that the bulk of star formation at z < 1 occurs in relatively mas- 
sive galaxies and in halos of mass 10" 5-12 5 Mq. 



results imply that a typical star in the Universe forms in galax- 
ies of similar mass, both in terms of stars and dark matter, 
from z ~ 1 to z ~ 0. The principle difference with redshift 
is that the distribution of masses increases with time, so that 
stars are more likely to form in a variety of systems at later 
epochs. 

Finally, Figure [TT] compares our model predictions for the 
evolution of the cosmic SFR density and stella r mass den- 
sity to data compil ations provided by iHopkinsI d2004l) and 
IWilkins et alj (H008), respectively. The top panel includes the 
model prediction fo r the SFR under the two different assump- 
tions discussed in $2.6\ The agreement between model and 
data at z < 1 in this figure is largely by construction, as men- 
tioned in £12.31 but is included here for completeness. The 
discrepancy between model and data at z > 1, which is more 
apparent in Figure Q] but is also seen here, is discussed in the 
next section. 

3.5. The model at z> 1 

In this section we have focused largely on redshifts less than 
one. At higher redshifts the model fails to match the observed 
cosmic SFR stellar mass density evolution at z > 1 as shown in 
FigureQj] and, relatedly, the normalization of the stellar mass 
function at z = 2, shown in FigureQ] This disagreement arises 
due to a more generic discrepancy between observed SFR in- 
dicators and stellar mass estimates, as discussed in $2.3\ As 
discussed in that section, one possible explanation is that the 



FIG. 1 1 . — Top Panel: Cosmic SFR density as a function of redshift. Our 
model predictions both for the no-merger (solid li ne) and merger (dashed 
line) models are compared to the data compilation of Hopkins 12004, circles). 
Bottom Panel: Cosmic stellar mass densit y as a function of redshift (solid 
line) compared to the data compilation of Wilkins et al. 12008, circles). 



IMF evolves with redshift (e.g. lWilkins et alj|2008l) . although 
we emphasize that this possibility is controversial. Nonethe- 
less, the generic inconsistency between observations implies 
that our model cannot be robustly applied to z > 1 . 

Moreover, it is plausible that our assumption of a tight cor- 
relation between stellar mass and halo mass breaks d own at 
higher redshift (cf. discussion in lConroy et al.ll2008l) . This 
tight correlation is strongly supported at < z < 1 by the 
observed stellar mass-dependent autocorrelation function of 
galaxies, in the se nse that more massive galaxies are more 
strongly clustered (iLi et alJl2006l:lMeneux et al.l l2008). This 
observational result can be most easily explained if more mas- 
sive galaxies reside in more massive halos because halo clus- 
tering stren gth is a monotonically increasing func tion of halo 
mass (e.g. lZehavi et aU2005l : IConrov et al. 2006). 

This observat ional trend has not bee n unambiguously con- 
firmed at z > 2 (Adelberger et al. 2005), except perhaps at the 
very highest masses (Ouadri et al. 2007). It is clear however 
that there is strong restfram e UV luminosity-dependent clus- 
tering at these early epoc hs (I Adelberger et al.f2005HLee et alJ 
2006; Ouc hi et afll2005l) . By analogy with stellar masses at 
low redshift, this trend can be understood if UV luminosity, 
and hence the star formation rate, is strongly and monotoni- 
cally correlated with dark matter halo mass. If this is the cor- 
rect interpretation, th en our model must be modified at these 
early epochs (see e.g. lConrov et al.ll2008l) . 

3.6. Dependence on cosmological parameters 

The halo mass functions and halo mass accretion histories 
in this model are dependent on cosmological parameters. Be- 
cause it affects the shape and normalization of the mass func- 
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tion, the normalization of the power spectrum can have a large 
affect on our results. T he a nalytic framework for these halo 
properties described in £j2.2l allows us to straightforwardly ex- 
plore the effect of cosmological parameters on our results. 
Here we just consider the effect of the normalization of the 
power spectrum as specified by the rms fluctuations measured 
in 8 Mpc/i -1 spheres, erg. 

We find that the impact of erg on our results is impercep- 
tible for z = halo masses less than ~ 1O 14 M . This is due 
to the fact that the dependence of the halo mass function on 
erg is much stronger at the massive end. When considering 
a change from our fiducial model with erg = 0.76 to a model 
with erg = 0.90, even at M vir = 10 14 M Q the difference in halo 
abundance is < 0.3 dex, and at M vn - = 10 13 M Q it is < 0.1 dex. 
The accretion histories are also a function of erg, but again the 
effect is only manifest at high halo masses. Since the bulk of 
our results focus on halo masses < 10 14 M©, we conclude that 
the uncertainty in erg does not impact our conclusions. 

4. DISCUSSION 
4.1. Downsizing: what, when, and where 

The phenomeno n known as "downsizing", coined by 
ICowie etaf] (119961) . has received much attention recently, 
and, perhaps confusingly, has been attributed to a number 
of related but nevertheless different phenomena. In its most 
general sense the term highlights a shift in a preferred mass 
scale of a phenomenon related to stellar growth or star for- 
mation. With an observationally-constrained model for the 
redshift-dependent connections between star formation, stel- 
lar mass, and halo mass, we are in a position to clarify and 
outline the relations between the various meanings of down- 
sizing. For clarity, we focus discussion on galaxies with 
M s tar < 10 u Mq and z < 1, wher e our results are most reli- 
able (see e.g. lCattaneo et a l. 2008, for a theoretical discussion 
of downsizing for higher mass systems). 

Originally, downsizing described the observation that the 
maximum /T-band luminosit y of galaxies abo ve a SSFR 
threshold decreases with time C Cowie et al.ll9 96). In this def- 
inition, the SSFR threshold is independent of redshift. In Fig- 
ure|7]we can see that a line of constant SSFR will indeed in- 
clude more massive galaxies at earlier epochs, consistent with 
this notion. However, inspection of the full relations in this 
figure shows that this notion of downsizing is driven by the 
global phenomenon that all galaxies have lower star forma- 
tion rates at later times. In fact, these relations do not appear 
to show any preferred scale with stellar mass (except possibly 
at very high stellar masses), but rather they shift self-similarly 
in time as noted by Noes ke et al.l (l2007blah . 

Figures [8] and [9] display another type of downsizing in the 
sense that the dark matter halo mass at which star formation is 
most intense shifts to lower masses at later times. This trend is 
apparent in both our favored model where stellar growth is en- 
tirely due to star formation (solid lines) and in the model that 
includes stellar growth due to mergers (dashed lines), sug- 
gesting that this form of d ownsizing is a generi c feature of 
dark matter halos. Indeed, iNeistein et al.l (120061) has argued 
that downsizing arises naturally from the accretion histories of 
the dark matter halos themselves. While this is an intriguing 
possibility, the uncertain relation between halos and galax- 
ies (connected in their terms by the competition between gas 
heating and cooling), makes their conclusions difficult to in- 
terpret at face value. 

Downsizing has also been attributed to the observation that 
more massive galaxies seem to have formed the bulk of their 



stars earlier. This type of downsizing has been referred to 
as "archaeological downsizing" because it is observed in the 
fossil record of the spectra of z ~ galaxies. It has been 
most convincingly demonstrated in local elliptical galaxies 
where one finds that more massive galaxies formed the bulk 
of their stars e arlier and over short er timescales than less mas- 
sive galaxies (Thomas et al. 2 005b . 

This form of downsizing can be seen clearly in Figure [4] 
which shows that the most massive galaxies formed the bulk 
of their stars earlier than less massive galaxies. From Figure|6] 
it is also clear that the peak of the SFR occurs at earlier times 
for more massive systems. Note however that in the top panel 
of Figure|6]there is no clear scale in the SSFR except perhaps 
for the most massive galaxies, and thus there is little evidence 
for any type of downsizing in this relation. 

A final meaning of downsizing concerns the sites where the 
bulk of stars are being formed at any epoch. This form of 
downsizing implies that stars are being formed in preferen- 
tially smaller systems at later times. In Figure [10] it is clear 
however that the typical masses hosting the bulk of star for- 
mation has not changed appreciably since z ~ 1 . There is thus 
no evidence for this form of downsizing given the available 
data. 

It is understandable, in light of the preceding discussion, 
that the term downsizing has been used to describe so many 
related but different phenomena, and that some authors find 
no evidence for a downsizing phenomenon. As we have seen 
in the various relations between SFR, SSFR, stellar mass, halo 
mass, and redshift, some show a shift in preferred scales with 
time, and some do not. Downsizing, of whatever type, thus 
manifests itself only in certain relations, and not in others. 

4.2. A characteristic halo mass? 

Recently several theoretical studies have raised the pos- 
sibility of a characteristic halo mass below which star for- 
mation occurs, and above which star forma t ion is truncated 
(iKeres et al.l 120051: iDekel & B irnboiml 120061: iBirnboim etail 
120071: ICattaneo et al.1 120071) . In this section we focus on the 
observational evidence for or against a sharp or narrow range 
in halo masses over which galaxy properties, such as star for- 
mation rates, change dramatically. 

This characteristic mass scale, which is thought to be ~ 
10 12 M© at z ~ 0, may be related to the observed stellar 
mass scale at which many properties of galaxies qualitatively 
change (e.g. iKauffmann et al.1 [20031) . The fundamental gas 
dynamical effect occurring in halos above this characteristic 
scale is thought to be the formation of a stable shock through 
which infalling gas must cross, thereby raising the tempera- 
ture of this newly accre ted gas to the virial temperature of the 
halo dKeres et al.ll2.005b . Accreted gas that is shock-heated is 
known as 'hot-mode' accretion, while gas that is not shock- 
heated is referred to as 'cold-mode' accretion. 

There are a number of outstanding issues related to any pos- 
sible sharp transition in e.g. galaxy colors or star formation 
rates occurring at a characteristic halo mass. First, the tran- 
sition from cold to hot-mode accretion does not a ppear to be 
partic ularly sharp in hydrodynamic simulations (Keres et al. 
2005). There is clearly a transition region, but it is broad, 
spanning the mass range ~ 10 1i_12 Mq. Moreover, the estab- 
lishment of a hot atmosphere does not guarantee that star for- 
mation will cease because such gas will still radiate and can 
thus cool (although hot, low density gas is more susceptible 
to further heating processes than cool, dense gas). Indeed, 
the cooling time of the intracluster medium at the centers of 
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massi ve clusters is in many cases < 10 9 Gyr dSanderson et al.l 
2006). One thus must propose additional mechanisms that 
are capable of supplying sufficient energy to keep the hot 
atmosphere from cooling and hence forming stars. Possi- 
ble mec hanisms include feedback from active galactic nu- 
clei (e.g. lCroton et al . 2006; Cattan eoet al.ll2008|). heating bv 
dynamical friction dMillerill986l:lKhochfar & Ostrikerl l2008l). 
thermal cond uction (Zakamska & Naravan 2003), virializa- 
tion heating ([W ang & Abel 2008), heating by ram-pressure 
drag (iDekel & Birnboim 2008), and supernovae heating. The 
relevance of these or other mechanisms to the shut-down of 
star formation is curren tly a subject of active debate (see 
Conroy & Ostriker 2008, for a recent evaluation). 

Our results on the relation between star formation rates and 
halo masses can shed light on this issue. In particular, our 
model provides a bridge between the observations and the un- 
derlying dark matter structure. At stellar masses < 1O U M0, 
where our results are most reliable, we find no significant ev- 
idence for a sharp characteristic halo mass at which star for- 
mation rates dramatically change, when considering average 
relations between star formation, stellar and halo mass. This 
statement is based on the following inferences. 

The scale at which galaxy properties such as color and mor- 
phology appear to chang e qualitatively is at M star ~ 10 10 3 M 
dKauffmann et al.l 120031) . Our abundance matching results 
shown in Figure [2] demonstrate that this stellar mass corre- 
sponds to a halo mass of ~ 10 I2 M Q at z ~ 0, in qualitative 
agree ment with the characteris tic halo mass scale mentioned 
above (Dekel & Birnboirnll2006l) . 

In Figure [8] however, for our favored model (the no merg- 
ing model) there is no abrupt change in the average star for- 
mation rate as a function of halo mass for M vir < 1O 13 M , 
at either z ~ or at higher redshifts. Instead, over this range 
the SFR-M v i, relation is approximately Gaussian with a broad 
peak at M vir < 10 12 ' 5 M Q at z ~ 0. We reiterate that the z — 
SFR-M vir relation is determined by 1) the z ~ SFR-M stal 
relation, where the model and data agree well, and 2) our con- 
nection between stellar and halo mass at z ~ 0. This latter con- 
nection is know n to reproduce the observed clustering p roper- 
ties of galaxies (IConrov et alj|2006t IZheng et~aT1l2007bl) and 
also agrees with weak lensing measurements of halo m asses 
as a function of stellar mass (Mandelb aum et af]|2006l) . At 
higher redshifts the peak shifts to higher halo masses, and the 
gradual roll-over at the high-mass end seen in the z ~ aver- 
age SFR-M v i r relation disappears. The data is thus consistent 
with there being no drop in star formation whatsoever above 
a given halo mass scale at higher redshifts, at least for halos 
with mass Mvb- < 1O 13 M , where we focus our results. In this 
figure the merger model does indeed produce a sharp break 
in the SFR-M v i r relation, but recall that this model fails to 
reproduce the observed SFR-M stal relation shown in Figure 

m 

It is important to stress that these conclusions at z ~ 
rest on the reliability of the o bserved z ~ a verage 
SFR -M star relation, as reporte d by ISalim et ail d2007l) (see 
also ISchiminovich et al J 12007). These authors caution that 
the star formation rates inferred for massive galaxies, M stal > 
10 11,5 Mq, may in some cases be upper limits because low 
levels of UV flux may arise from old stellar populations (e.g. 
iRich et all [2005) that are not typically included in modeling 
of star formation rates. Interpreting low levels of UV flux has 
historically been challenging for this reason. Similar issues 
arise at higher redshifts. Note however that we do not rely 
on these massive galaxies for our conclusions because they 



reside in very massive, M v [ r > 10 14 M Q , halos. 

Moreover, as shown in Figure [8] there is clearly no scale in 
the specific star formation rate as a function of halo mass — 
it is approximately a power law that scales as SSFRoc M~^ 5 
over at least two orders of magnitude in halo mass. Again, 
these statements apply to halos with mass M vr ,- < 10 13 M©. At 
higher masses our model is not well-calibrated. 

In sum, while a well-defined characteristic halo mass, above 
which star formation is truncated, may be an a ppealing mech- 
anism for generating red sequence galaxies (Cattan eo et al.l 
2008), there is no clear indication from our data-driven model 
that this scale is particularly sharp. It is clear that observed 
galaxy properties change qualitatively around a stellar mass 
scale of M star ~ 10 Mq, corresponding in our model to a 
halo mass of ~ 10 12 M Q . We simply emphasize that the data 
favors a rather gradual shift in galaxy properties across this 
halo mass scale. 

4.3. The relative importance of star formation and merging 
to galactic growth 

In order to translate our model predictions for stellar mass 
growth into predictions for star formation rates, we have to 
make assumptions for the fr actio n of mass growth attributed 
to mergers, as described in £12.61 One approach is to assume 
that all stellar material accreted onto the halo remains in the 
halo as satellite galaxies or is stripped and remains in the stel- 
lar halo. This is the no-merger model described above. In this 
model all stellar growth is due to star formation. The second 
approach is to assume that all of the accreted material imme- 
diately falls onto the central galaxy and hence contributes to 
its stellar growth. These two approaches should bracket the 
range of possibilities, as in reality some accreted material will 
lose energy and merge with the central galaxy, while other 
material will remain as bound satellites, or will merge with 
the central galaxy yet be dispersed outside the photometric 
radius. 

For galaxies with stellar mass < 1O IO M0, these two treat- 
ments for the importance of merging on stellar growth lead 
to indistinguishable predictions for the resulting star forma- 
tion rates (see e.g. Figure |7). Thus, we can state with con- 
fidence that galaxies below this mass range grow almost en- 
tirely by star formation, at least since z < 1 where we focus 
our analysis. This result can be understood as follows. Halos 
grow via the accretion of smaller halos. By inspection of the 
lower panel of Figure |2] it is clear that for halos with mass 
< 10 1 Mq, corresponding to stellar masses < 1O 1O M0, the 
fraction of available baryons that have been converted into 
stars drops precipitously. In other words, for these low mass 
halos, the even smaller mass halos that are contributing to halo 
growth are almost entirely devoid of stars. Furthermore, low 
mass halos h ave l argely completed their growth by z ~ 1, as 
discussed in £12.51 and thus any resulting stellar growth since 
Z ~ 1 must come from within the hal o, i.e. via star form ation. 
These points were also discussed in iPurcell et al.l d2007l) and 
are qualitatively consist ent with current semi-analytic mod- 
els dGuo & Whitdl2008l) . They robustly follow from the inte- 
grated star formation efficiencies shown in Figure [2] 

At larger stellar masses the two treatments yield different 
predictions for the star formation rates of galaxies. The re- 
sults presented in £13.31 show that the data on the SFR-M star 
relation match the approach that attributes all stellar growth 
to star formation, at least for stellar masses < 1O U M0 and 
z < 1, where we focus our analysis. Consideration of the mass 
range 10" < M star < 10 1L5 M Q , where our results are some- 
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what less certain (for the reasons discussed in 81,3.21 ). provides 
further evidence in favor of this 'no-merger' scenario, i.e., that 
incoming halo mergers do not contribute substantially to the 
growth of the central galaxy. It thus appears that over this 
entire stellar mass and redshift range, stellar mass growth in 
galaxies is dominated by star formation. These conclusions 
are largely consistent with results from cosmological hydro- 
dynamic simulations and may help explain the dominance of 
disk galaxies at these stellar masses, if disks are a signpost of 
a relatively quiescent history (Mailer et al. 2006). 

At first glance this may seem surprising because at these 
higher masses one expects accretion of halos massive enough 
to host large galaxies. The accretion of such objects is, as 
mentioned above, a generic prediction of A^-body simulations 
coupled to our connection between galaxies and halos. Of 
course, the accretion of stellar material onto the halo need not 
necessarily lead to growth of the galaxy residing at the center 
of the halo because the accreted material may either remain 
in orbit within the halo or may be tidally disrupted before 
it can spiral into the center. In the latter case, the material 
will contribute to the observed diffuse i ntracluster light that is 
ubiquitous in large dark matter halos dGonzalez et al.ll2005t 
Zibe ttiet all 12005). Indeed, our results indicate that some 
combination of these two sce narios is precisely w hat is hap- 
pening (see also discussion in lConrov et al.l20 07b). Evidence 
for the former scenario, whereby accreted material remains 
as bound satellites, is corroborated by the observed increase 
since z ~ 1 in the fraction of galax ies at a given halo mass that 
are satellites (IZhenget al.ll2007bl) . 

In sum, our results suggest that stellar growth since z ~ 1 
is dominated by star formation, as opposed to mergers, for 
stellar ma sses < lO^'M^. This conclusion echoes the conclu- 
sions of iBell et al.1 (|2007) who used the observed, redshift- 
dependent, star formation rate — stellar mass relations to 
'predict' the evolution of the stellar mass function since z = 1 . 
These authors then compared this predicted evolution to the 
actual evolution of the mass function in order to conclude that 
mergers had a minor effect on the growth of intermediate and 
low mass galaxies. 

5. SUMMARY 

This paper presents a model for the evolution of galax- 
ies that is based on the observationally-motivated assumption 
of a tight correlation between galaxy stellar mass and dark 
matter halo mass. This assumption is used to populate ha- 
los with galaxies from z = 2 to z = using theoretical halo 
mass functions and observationally-constrained galaxy stel- 
lar mass functions. Halos (and the galaxies within them) are 
evolved forward in time using estimates for halo growth cal- 
ibrated against A^-body simulations. This then provides the 
average stellar mass growth of galaxies as a function of z = 
stellar and halo mass. At M star < 10 10 Mq the model robustly 
predicts that the vast majority of stellar growth is due to in 
situ star formation since small halos do not accrete significant 
amounts of stellar material. At higher masses, where the halo 
merger rate is higher, mergers and accretion could in principle 
contribute to stellar growth. However, the model agrees with 
an array of data when all stellar growth at these higher masses 
is attributed to star formation (rather than a substantial frac- 
tion being due to mergers) for galaxies with M star < 10 u Mq 
at z < l • We do not use our model to address the growth his- 
tory of more massive galaxies because various aspects of the 
model become uncertain in this regime. 

With the assumption of a one-to-one correlation between 



stellar and halo mass, the only freedom within our framework 
is the particular form adopted for the redshift-dependent stel- 
lar mass function. We have adopted a form that provides the 
best fit to a variety of data including the observed stellar mass 
function at < z < 1, the cosmic SFR and stellar mass den- 
sity, and the SFR-M stal relation over the range < z < 1 . This 
model can thus be thought of, in part, as a self-consistent syn- 
thesis of the available data relating galaxy SFRs and stellar 
masses across time; it allows us to connect galaxy populations 
at a given epoch with those at another epoch. The model also 
effectively connects the observations to the underlying dark 
matter structure, thereby providing a bridge between obser- 
vational results and theoretical work aimed at understanding 
such observations. 

The principle new result that can be obtained from this 
framework is a directly-constrained form of the star forma- 
tion rate in galaxies as a function of halo mass. Our approach 
provides a direct link between observations and these mod- 
els in the sense that any model which reproduces this con- 
strained relation for SFR(M V ;, -,z) will automatically match the 
wide variety of observational results discussed herein, over 
the last half of the Universe's age. This result can thus help to 
distinguish between the processes responsible for triggering 
and halting star formation, and can be directly employed in 
constraining models and simulations of the physics of galaxy 
formation. 

The success of this simple model at describing an array of 
data over the stellar mass range 10 9 < M star < 10 n M Q and 
redshift range < z < 1 indicates that the relation between 
galaxies and halos is surprisingly simple, smooth, and mono- 
tonic over these ranges. The most significant short-coming of 
this model is its inability (in its present form) to predict dis- 
tributions of properties, rather than averages, as a function of 
stellar and halo mass. Such information is clearly needed to 
understand the color bi-modality seen in the color-magnitude 
diagram, as well as the detailed properties of satellite galax- 
ies, and we will address this in future work. 

This model relies on observational inputs that are rather un- 
certain, such as the evolution of the stellar mass function and 
the IMF, and the quantitative predictions of this model are 
thus necessarily uncertain. Despite these unavoidable uncer- 
tainties, the general trends predicted by this model, such as the 
dependence of the SFR of galaxies on galaxy and halo mass, 
are robust and highlight the underlying connections both be- 
tween the panoply of observations at high and low redshift, 
and between the observations as a whole and the underlying 
dark matter distribution. 
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